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ABSTRACT 

We present a catalogue of molecular cloud cores drawn from high latitude, medium 
opacity clouds, using the all-sky IRAS Sky Survey Atlas (ISSA) images at 60 and 
100 jim.. The typical column densities of the cores are N(H2) ^ 3.8 x lO^^cm"^ and 
the typical volume densities are n(H2) 2 x lO^cm"'^. They are therefore significantly 
less dense than many other samples obtained in other ways. Those cloud cores with 
lists' point sources are seen to be already forming stars, but this is found to be only a 
small fraction of the total number of cores. The fraction of the cores in the protostellar 
stage is used to estimate the prestellar timescale - the time until the formation of a 
hydrostatically supported protostellar object. We argue, on the basis of a comparison 
with other samples, that a trend exists for the prestellar lifetime of a cloud core to 
decrease with the mean column density and number density of the core. We compare 
this with model predictions and show that the data are consistent with star formation 
regulated by the ionisation fraction. 

Key words: stars: formation ~ ISM: globules. 



1 INTRODUCTION 



The denser regions of the interstellar medium are usually 
known as molecular clouds, and stars form in the very high- 
est density parts of these clouds, which are usually referred 
to as molecular cloud cores. The process by which material 
is turned from a molecular cloud core into a star is far from 
being fully understood, since many complex physical pro- 
cesses are involved (see, for example, Mouschovias 1991, for 
a review). 

The key to understanding the process of star forma- 
tion appears to lie in the earliest stages of molecular cloud 
core evolution, since the initial conditions determine much of 
what subsequently takes place. Many observational studies 
have been conducted of such dense cores (e.g. Myers & Ben- 
son 1983; Myers, Linke & Benson 1983; Clemens & Barvainis 
1988; Benson & Myers 1989; Bourke, Hyland & Robinson 
1995a), usually starting from an optically selected sample of 
dark clouds on sky survey plates. In this paper we present 
a sample of cloud cores selected on the basis of far-infrared 
optical depth using the IRAS Sky Survey Atlas (ISSA) im- 
ages (Wheelock et al. 1994), with a view to broadening the 
range of the physical parameters of cores that have been ex- 
plored, and hence gaining further insight into the evolution 
of molecular cloud cores. 

The paper is laid out as follows: Section 2 describes 
the manner in which optical depth maps were constructed 
from the ISSA data, including technical details such as back- 
ground subtraction, and goes on to discuss cloud core se- 



lection and methodology verification techniques; Section 3 
describes the new molecular cloud core catalogue that we 
have constructed from the data and the associations with 
previously known molecular clouds and IRAS point sources; 
Section 4 contrasts the properties of the new catalogue with 
previous molecular cloud core catalogues and compares the 
ensemble of mean catalogue properties with the theoretical 
predictions of models of ionisation-regulated star formation; 
and Section 5 presents the conclusions of the paper. Read- 
ers who are more theoretically inclined may wish to read 
Sections 3 & 4 before going back to Section 2. 



2 OBSERVATIONAL DETAILS 
2.1 Background subtraction 

We chose to select the clouds using the IRAS Sky Survey 
Atlas (ISSA), which is an all-sky set of images at each of 
the four IRAS wavebands of 12, 25, 60 & 100 fim (Beich- 
man et al. 1988; Wheelock et al. 1994). However, since we 
wished to select the coldest, densest clouds, we concentrated 
on the two longest wavelengths, 60 & 100 /xm. We used these 
to construct optical depth and temperature maps of the re- 
gions of interest, in a similar manner to that adopted by 
Wood, Myers & Daugherty (1994). These latter authors se- 
lected previously known molecular cloud regions, whereas 
we endeavoured to construct a distinct sample of previously 
little-studied clouds drawn from the all-sky set of plates. But 
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before optical depth and temperature maps can be made, it 
is necessary to ensure that all background emission has been 
removed from the images. This is because in making opti- 
cal depth and temperature maps, one must take the ratio 
of emission at two different wavelengths (in this case 60 & 
100 /xm). Therefore, any offset at one or other wavelength 
due to background emission not associated with the object 
will affect the ratio measurement. The ISSA images have al- 
ready been fully processed and corrected for most of the Zo- 
diacal background emission, but they still contain extended 
emission from the Galactic Plane and some residual Zodiacal 
emission. These are the backgrounds we now address. 

Pixel histograms of all the ISSA fields were constructed 
by binning the pixels in each field into histograms of sur- 
face intensity. On examination it was found that some fields' 
pixel histograms contained a single peak, some contained 
double or multiple peaks and some had much more compli- 
cated structures. Each field was searched to identify whether 
the histogram had a peak with a width of less than 2MJysr~^ 
(i.e. ~ 10 X the calibration noise) and at low enough inten- 
sity to be consistent with an area of background. Any such 
peak was taken as evidence of an area of low level emission 
in the field which can be described as an area of background 
containing only low level cirrus. A number of the fields were 
inspected visually. It was found that the pixels in the low 
level peak of the histogram in each case came from a con- 
tiguous, discrete area of the sky, and were not randomly 
isolated pixels, confirming that this automatic method was 
indeed finding genuine areas of sky. 

A similar method was previously carried out for opti- 
cal images by Beard, MacGillivray & Thanisch (1990), who 
showed that histograms of pixels' intensities contain valu- 
able information about the background regions within the 
image. In relatively empty areas of the sky the pixel his- 
togram takes the form of a single Gaussian peak whereas in 
regions densely populated with real sources several further 
peaks appear in the histogram at higher intensities. 

The position of the peak was recorded and used as a 
measure of the background surface brightness in the field. 
272 ISSA fields were identified as having suitable back- 
ground regions by this technique at 100/im and 368 were 
found at 60/xm from a total of 430 for each wavelength. The 
widths of the low intensity peaks in the pixel histograms 
were generally at least a factor of two or three larger than 
the average calibration errors for the field - as estimated 
from the standard deviation in pixel values from one Hours 
Confirmed scan (HCON) to the next - indicating the exis- 
tence of real cirrus structure within the background regions. 

To investigate the overall nature of the background, 
maps which interpolated from region to region of back- 
ground were constructed. This was done by first recording 
the values of the pixel histogram low intensity peaks (in 
MJysr"^) and recording the central position of the field 
in which it was found. A simple interpolation technique 
was devised, somewhat akin to box-car averaging. For ev- 
ery square degree on the celestial sphere the distance to the 
centre of the nearest few ISSA fields containing a region of 
background was calculated, to allow averaging. This gives 
a series of angular displacements to the background regions 
(^1, 02,^3, 6*4, ....^n), where: 

cos(Si) = cos(deci)cos(dec)cos(ra — rai) — sin(deci)sin(dec). 




Galactic 



Figure 1. Interpolated map of 100-fim background regions. Val- 
ues range from —0.2 to MJySr~^. There is a gap in the data 
through the map along the Galactic Plane (partly obscured be- 
cause of the projection), due to a larger number density of clouds 
near the Galactic plane causing source confusion and thus pre- 
venting the background subtraction method from working. 



Here the right ascensions and declinations of the positions to 
which we wish to interpolate are denoted ra and dec respec- 
tively, and the position of each background region is denoted 
rai and deci. We then attributed to each position the value 
of surface brightness: 

E^-"^cos(^)/. 
^^■<^cos(^) • 

This was chosen because it interpolates between positions 
smoothly, does not produce discontinuities and weights more 
heavily the nearest background regions, even given the non- 
uniform spatial sampling of background surface brightness. 
It essentially involves smoothing by a cosine bell of radius 
0. The result of this process was then projected onto a 2d 
surface with equal area projection to create a map of back- 
ground brightness for the whole sky. 

Fig. 1 shows the map of the background constructed by 
this technique at lOO^m for (j) — 12 degrees. The plot is an 
'equal area' projection in Galactic coordinates. The intensity 
varies smoothly, and increases towards the Galactic Plane in- 
dicating that the cirrus intensity within background regions 
increases towards the Plane. The map contains negative val- 
ues showing that the Zodiacal background subtraction used 
in producing the ISSA dataset was liable to over-compensate 
in some regions. The apparent gap through the centre of the 
map is due to our having to discard regions of high source 
confusion in the Galactic Plane, where our Galactic back- 
ground subtraction did not work due to source confusion. 

The residual 'striping' characteristic of images produced 
by IRAS was reduced to very low levels in the ISSA images, 
because of the careful calibration used to produce them. 
However, the act of ratio-ing two images tends to enhance 
any striping effects that remain. In some fields this striping 
was clearly visible, and in some cases dominated the struc- 
ture of the optical depths maps. These fields were discarded. 
We selected 60 of the fields used in the interpolation for fur- 
ther study, based on their background regions being clearly 
identifiable and measurable, and on the quality of the im- 
ages. 
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Table 1. List of cloud positions, associations, velocities and distances. In column 4 a superscript indicates that the association is a: (1) 
dark nebula catalogued by Lynds (1962); (2) bright nebula catalogued by Lynds (1965); (3) molecular cloud identified and mapped in 
CO by Taylor, Dickman & Scovillo (1987); (4) cluster listed by Lang (1992); and (5) molecular cloud identified and mapped in CO by 
Ramesh (1994). 



Cloud 


R.A. 


Dec. 


Associations 


Velocity 


Distance 


Name 


(1950) 


(1950) 


identified 


(LSR)/kms-i 


(pc) 


021A 


OO'' 


00"' 


— 76° 


00' 


Chameleon 


- 


200 


270A 


OO'' 


02"' 


+13° 


50 


None 


- 


<80 


422A 


Ol'' 


30'" 


+76° 


30' 


Cepheus 


- 


300-800 


423A 


02*^ 


10'" 


+75° 


30' 


Cepheus, L133\ LBN^, 


+3.1^ +3.5 


800 












T4863, T486'*, Berkeley8* 






002B 


02'' 


30'" 


-85° 


00' 


Chameleon 


- 


200 


423B 


03'' 


00'" 


+81° 


20' 


Cepheus 


- 


300-800 


411A 


04^ 


15™ 


+64° 


20' 


Cepheus 




ouu ouu 


205A 


04^ 


28™ 


+04° 


00' 


Orion 




450 


050A 


oyh 


26™ 


-48° 


00' 


R42^, Merlotte 66* 


+4.8 


<250 


221A 


15'' 


44™ 


-04° 


00' 


p Ophiuchus, L134N1 


+2.6 


125 


420C 


20'' 


40™ 


+67° 


00' 


Cepheus, L11481, T3793 


+3.5 


300 


420B 


21'' 


30™ 


+66° 


00' 


Cepheus, L11761, T3093, 


-10.8 


800 












NGC 7142"*, NGC 7129* 






002A 


21'' 


30™ 


-83° 


00' 


Chameleon 




200 


334A 


21'' 


48'" 


+35° 


30' 


None 




<270 


422C 


22'' 


00™ 


+77° 


30' 


Cepheus, LBN^, 


-1.7-+ +5.1 


300 












T4203,T4283,T43l3 






422B 


22'' 


45™ 


+73° 


40' 


Cepheus, T'i2&\ T43l3 


-3.7 


300 



2.2 Colour temperatures 

We used the STARLINK data reduction package IRAS90 
and, in particular the routine COLTEMP, to create colour 
temperature and optical depth maps from the ISSA images. 
We briefly describe the technique here (for further details, 
see Berry 1993a & b). 

The far infrared radiation from a cloud of temperature 
T emitting a black body spectrum, B(i/, T), and absorbing 
I(:^, T)dr^, at a position in the cloud with an optical depth 
Tij, leads to a flux received by an observer f{iy,T) given by: 

f(i.,T) = {l-e'^")B{u,T). 

Generally t^, is dependent on u in such a way (Hilde- 
brand 1983) that: 

where f3 is the dust emissivity index and i^c is the critical 
frequency at which the optical depth is unity. Throughout 
this work we use the value suggested by Hildebrand (1983) 
of /3 = 1 at far-infrared wavelengths. These two equations 
describe a 'greybody' spectrum. We assume that the cloud 
is optically thin at the wavelengths observed - Wood et al. 
(1994) show that even towards the Galactic Plane this is 
true. Because r,y <C 1, 1 — e""^" ~ and we therefore use 
the expression: 

f(^,T)~ (^)b(/.,T). 

The IRAS detectors were sensitive over a wide bandpass 
and there were 4 separate wavebands (i=l, 2, 3, 4 for 12, 
25, 60 and 100 /xm respectively) so that the measured flux 
in waveband i is: 



fi = y Ri{u) (^-^^ B{u,T)du 

where Ri(/^) is the spectral response curve for the waveband 
i receiver (see Beichman et al. 1988). 

By taking the ratio of intensities at two wavebands i 
and j, and using the preceding equation, one can derive 

fi _ /Ri(i^)!/B(z^,T)dj/ 
f^ " J Ri{u)uB{u,T)du 

This value is dependent on T, and the response curves of 
the receivers. Using the listed response curves in the IRAS 
Explanatory Supplement (Beichman et al. 1988), one can 
tabulate fi/fj versus T. In the routine COLTEMP a spline 
giving T as a function of fi/fj is created by fitting to the 
tabulated values. For any observations of a cloud at 2 wave- 
lengths one can then estimate the temperature and calculate 
the critical frequency at which r^^ = 1: 



Using this, one can calculate the optical depth of the 
cloud at another wavelength by using the equation for r 
above. We altered COLTEMP to allow temperatures as low 
as lOK to be used (the publicly available version only accepts 
temperatures greater than 30K). In this way, optical depth 
maps were made of our chosen regions at 100 /im. 



2.3 Optical depth maps 

The optical depth at 100/im due to dust along the line of 
sight is expressible (Hildebrand 1983) as: 

Tioo = 7r(a)^QiooA^g, 
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Figure 2. (a) 100-/^m optical deptli map of tlie molecular cloud L1689, produced from the ISSA sky survey plates. Base 
contour level is 0.8 xlO~^ and the contour interval is 0.3 xlO""'. (b) ^^CO map of the same field taken from Loren (1989). 



where A^g is the column density of grains, (a) is the average 
radius of the grains and Qioo is the emission efficiency of 
the dust grains at lOO/xm. The mass column density of the 
grains is then given by: 

where p is the average grain density. The total mass of the 
dust in a cloud of projected area A is therefore given by: 

Mdust = A{ad}, 

where (ad) is the mean column density of dust in the cloud. 
The typical dust-to-gas mass ratio in the local ISM is nor- 
mally assumed to be approximately 1:100. However it is clear 
that at 100 pm a significant fraction of the cold dust in 
the ISM is not detected, due to temperature and optical 
depth effects (see e.g. Wood et al. 1994). Hence a simple 
application of this ratio to these data will underestimate 
the total mass of gas. Wood et al. (1994) argued that only 
1/50 of the total dust is detected at IRAS wavelengths, and 
used a dust-to-gas ratio of 1:2000 in their analysis. They 
arrived at this value by using two relations: Firstly they 
adopted Av ~ 2 x 10'' x rioo (Langer et al. 1989); and 
secondly they used A^H2(cm~^) ~ 10"^^ x Ay (Bohlin, Sav- 
age & Drake 1978). We follow Wood et al. (1994) and use 
a value of 1/2000 in the current study. With a value for 
(ap/Qioo) ~ 32gcm~^, one obtains the expression for the 
mass of material in a cloud: 

M,,o„d/M0 ~ 1.25 X 10"' x I>'(pc) noo, 

cloud 

where D is the distance to the cloud, and the optical depth 
of each pixel in the cloud image is summed. 



2.4 Cloud selection 

Of the resultant optical depth maps, some contained from 
one to three cloud complexes, varying in size from a few pix- 
els to half a field, while several more were discarded, either 
due to having very little structure, being too noisy, or having 
major contamination from residual Zodiacal emission that 
none of the processing had been able to remove. Some of 
the clouds were found at the edges of the fields and hence 
were truncated. /.From the remaining maps, 17 clouds were 
selected for further study. The clouds we selected are listed 
in Table 1: Column 1 lists the name we assigned to each 
cloud, derived from the ISSA field number in which it was 
found; Columns 2 & 3 give the approximate position of the 
centre of the cloud; and the remainder of the table assigns 
distances, velocities and associations to the clouds we have 
selected. 

It was found that 3 of the clouds were previously iden- 
tified by Lynds as dark clouds (Lynds 1962), 2 contained 
Lynds bright nebulae (Lynds 1965), 7 had been identified 
by other authors (Taylor et al. 1987; Ramesh 1994), who 
had subsequently measured their CO velocities, and 3 had 
nearby open cluster associations. Of the 3 open clusters, 
two have been dated and were found to be old: NGC 7142 
is thought to be 4 billion years old; Merlotte 66 is 6 billion 
years old; in both cases implying that they were probably 
not linked with the cloud. In addition, comparison with the 
CO Galactic plane surveys (Dame et al. 1987) revealed that 
several clouds were associated with known cloud complexes. 
Only two had no previously published associations. Column 
4 of Table 1 lists the known cloud associations. 

Unlike molecular maps and surveys, which give the ve- 
locity of the clouds and hence give an estimate of the dis- 
tances, these ISSA selected clouds - like the optically se- 
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lected clouds of Lynds (1962) - do not have easily derivable 
distances. Distances were estimated either from velocity and 
spatial association with the Orion, Cepheus, Chameleon and 
Ophiuchus complexes, or by estimating an upper limit for 
distance obtained by assuming the clouds lie in the Galac- 
tic Disc - i.e. less than 60 pc away from the Galactic Plane 
(c.f. Clemens, Sanders & Scoville 1988). Clouds found in 
Cepheus present a particular problem when one attempts to 
assign a distance. There are two different complexes along 
the line of sight: one at approximately 300 pc with a veloc- 
ity kms~^; and one in the local spiral arm at approxi- 
mately 800 pc and with a velocity of ~ — 12kms~^ (Grenier 
et al. 1989). Some of the Cepheus clouds we selected had 
been sampled with CO observations (Taylor et al. 1987), 
and the measured velocities revealed the clouds belonged to 
one or other complex. This revealed that the sample pre- 
sented here contains clouds from both complexes and hence 
that the three Cepheus clouds in the sample without CO 
observations could be at either of the two distances. In to- 
tal, 11 of the 16 clouds had a single distance assigned to 
them, 2 had upper limits, and 3 could be at either of 2 dis- 
tances. Column 5 of Table 1 gives the velocity of the cloud 
(if previously measured), and column 6 gives the estimated 
distance. 

2.5 The case of L1689 

As a cross-check of both the procedure used and the mass 
estimates derived, we made a map of a previously studied 
cloud, L1689. Colour temperature and optical depth maps 
were constructed and the results are presented in Figure 2(a) 
as a 100-/im optical depth map of the region. An extended 
region of high 100-^m optical depth containing some struc- 
ture is seen. Figure 2(b) shows a ^''CO map of the same 
region taken from Loren (1989). A similar morphology is 
seen in both maps. The two star-forming cores R57 (alias 
L1689N), and R59 (alias L1689S) are clearly visible, as is 
the isolated core R65 (alias L1689B). 

At the centre of the cloud there is an apparent 'hole' 
in the optical depth map at the position of the bright point 
source, IRAS16288, which is a young protostar at the cen- 
tre of L1689S. Wood et al. (1994) also noted that a bright, 
point-like source dominating the IRAS emission can cause 
an apparent hole in the 100-/im optical depth. They ascribed 
this to a beam filling factor effect. The right hand side of the 
above equation for T(ioo) also contains a term for the solid 
angle of the source, fl, which has been set equal to 1 in our 
analysis. For an extended molecular cloud source this is a 
valid assumption, but it fails when there is a bright point 
source in the beam. This is the effect we see in the case of 
IRAS16288. 

We calculated the mass of L1689, using the above equa- 
tion, from the 100-pim optical depth map in Figure 2(a), 
and found a value of 448 Mq (assuming a distance of 160 
pc). Loren (1989) used the '^^CO data shown in Figure 2(b), 
and obtained a value of 566 Mq. These two measurements 
are consistent to within ±20 per cent, which is as accurate 
as either method can claim, and so we conclude that the 
technique outlined not only gives accurate qualitative infor- 
mation on the morphology of interstellar clouds, but also 
appears to provide relatively good quantitative estimates of 
cloud masses. Nonetheless, we realise that masses derived 



from 100-^m optical depths might be underestimated in 
some cases, if the 100-/xm emission is completely optically 
thick. In none of the cloud cores that we studied did this 
appear to be the case. 



3 THE CORE CATALOGUE 
3.1 Core properties 

Figures 3 & 4 show grey-scale images, with isophotal con- 
tours overlaid, of nine of our molecular cloud regions. They 
can be seen to vary in size from 1 to 5 degrees across, and also 
to vary in complexity and structure. Some residual striping 
can be seen in two of the fields. However the typical struc- 
tures associated with molecular clouds can be seen in all 
of the images - namely cores, filaments and other struc- 
tures on all scales within the maps. Some of the structure 
seen in the maps corresponds to previously known molecular 
clouds, but some of the clouds, and many of the cores, had 
not been previously catalogued. For example. Figure 3(d), 
420B, and Figure 4(d), 423A, both contain a Lynds dark 
nebula and an open cluster. 423A also contains an HII re- 
gion and two further catalogued clouds. Some of the regions 
contain no previous associations - especially those in the 
southern hemisphere. For a full list of known associations 
see Table 1 and section 2.4 above. 

A catalogue of the most opaque regions in each cloud 
was produced and the maps were examined. The definition 
of what constitutes a core in a molecular cloud is of necessity 
somewhat subjective. We chose to define a core as the most 
opaque 1% of a cloud's area. This level was chosen to ensure 
that the area defined as a core was relatively small compared 
to the cloud, and hence would only include the most dense 
regions where we would expect star formation to take place. 
We then calculated the mass of each core. 



Table 2 lists the properties of our new core catalogue. 
Column 1 lists the new number we have designated for each 
core in order of Right Ascension. Column 2 lists the cloud 
in which it is located, following the numbering convention 
of Table 1. Columns 3 & 4 list the position of the centre 
of each core. Column 5 lists the derived colour temperature 
of each core. Column 6 lists the solid angle of each core in 
square arcmin. Column 7 lists the mean optical depth, and 
column 8 lists the peak optical depth of the core. Column 9 
gives the ellipticity, column 10 lists the mean radius, column 
11 lists the position angle (north through east) of the major 
axis and column 12 lists the mass of each core. 

The median 100- /^m optical depth of our 60 cores is 1.9 
X lO"**, corresponding to a column density of N(H2) — 3.8 
X 10^^ cm~^. The median radius of the cores for which a 
distance is known is 0.31 pc, with a mean radius of 0.41 pc. 
The mean is skewed by a small number of large cores, so 
we prefer to use the median to characterise our sample. The 
median volume density of the sample is ~2 x 10^ cm~^ 
(the mean volume density is very similar). Hence we see 
that, by selecting our core sample based on a wavelength of 
100 fim, we have typically selected somewhat lower density 
cores than many previous surveys of molecular cloud cores 
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Figure 3. 100-Atm optical depth maps, with contours overlaid, of: (a) 002A; (b) 050A; (c) 205A; (d) 420B; (e) 420C. 
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Right Ascension (B1 950) 



20' 




W Right Ascension (B1950) 

Figure 4. 100-/im optical depth maps, with contours 
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(d) Rigi^^t Ascension (B1950) 

overlaid, of: (a) 422A; (b) 422C; (c) 422B; (d) 423A 
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Table 2. The catalogue of molecular cloud cores, together with their positions, sizes, optical depths, shapes, orientations and masses. 



JWT 


Cloud 




R.A. 




Dec. 




Temp 


Size 


T 


T-poak 


EUip- 


r 


e 


Mass 


No. 


name 


(1950) 


(1950) 


(K) 


(arcmin^) 


(xlO-3) 


(xlO-^) 


ticity 


(pc) 


n 


(Mo) 


1 


270A 


uu 


riom 
Uo 


Uo 


-rio 


OD 


OU 


29.2 


23.4 


0.02 


0.018 


0.43 


<0.2 


11 


<2 


2 


422A 


Ui 


ZL 


00 


+ /o 


/1 1 ' 
41 


( 


21.7 


105 


0.40 


0.439 


0.43 


0.5-1.2 


67 


300-2100 


3 


422A 


Ui 




04 


H- (u 


A a' 
44 


Z 1 


21.5 


45.3 


0.38 


0.404 


0.71 


0.3-0.9 


68 


120-880 


4 


423A 


Uz 


Uo 


io 


1 vco 
+ (0 


01 


OU 


20.8 


54.7 


0.52 


0.580 


0.68 


0.38 


71 


200 


5 


423A 




nnm 
ZU 


40 


-f- ( 


99/ 
ZZ 


04 


21.4 


12.5 


0.49 


0.505 


0.18 


0.18 


134 


44 


6 


002B 


uz 


01 m 
oi 


CCS 

DO 


QA 
— o4 


OD 


Uo 


23.8 


21.9 


0.12 


0.142 


0.01 


0.16 


123 


8.3 


7 


423A 


uz 


oD 


zU 


+ ( 4 


( 


Uo 


20.2 


26.6 


0.55 


0.651 


0.41 


0.26 


116 


100 


8 


423B 


noh 


OZ 


/I OS 
4z 


+ OI 


UU 


zU 


22.0 


17.2 


0.19 


0.206 


0.69 


0.2-0.5 


35 


20-160 


9 


423B 


Uo 


1 c;rn 

io 


CCS 

00 


H-oi 


OU 


10 


19.8 


54.7 


0.24 


0.411 


0.36 
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Table 3. IRAS point sources associated with the cores in this catalogue, together with their flux densities and other known associations. 
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(see below) . The difFerence is probably mainly due to our se- 
lecting relatively isolated clouds as a result of our constraints 
on background emission. 



3.2 PSC associations 

All of the IRAS point sources within the boundaries of each 
core were located using the Point Source Catalogue (PSC), 
and are listed in Table 3. The core in which the source is 
found is noted along with the source flux densities at each 
IRAS waveband. A superscript in the last column indicates 
that the IRAS PSC listed a known association for the source: 
(1) indicates that the source appears in the Smithsonian As- 
trophysical Observatory Star Catalogue (SAO); (2) repre- 
sents the catalogue of Ohio State University Radio Sources; 
(3) is the IRAS Serendipitous Survey Catalogue; (4) is the 
Dearborn Observatory Catalog of Faint Red Stars; and (5) 
is the IRAS Small Scale Structure Catalog (Beichman et al. 
1988 and references therein). 

There are 21 source associations, of which: 11 are de- 
tected only at 100 fim; two are detected only at 60 /im; with 
one detected at 60 & 100 fim, but no shorter wavelengths; 
five are only detected at 12^m; one source is detected at 
12 & 25 fim, but not at longer wavelengths; and one source 
is detected at all four wavebands. The eleven 100-/.im-only 
sources have upper limits at the shorter wavelengths such 
that we can say that their spectral energy distributions peak 
at a wavelength of around 100 /xm or longer. The one source 
detected at all four wavebands and the source detected at 
60 & 100 fim only, also have rising spectra to longer wave- 
lengths. The other two 60-/im-only sources are also consis- 
tent with having rising spectra. This is the form of spec- 
tral energy density we would expect for deeply embedded 
young stellar objects (YSOs) or protostars. The six sources 
detected only at the shortest wavelengths have spectra con- 



sistent with field stars or other objects. Hence there appear 
to be 15 PSC sources associated with our 60 cores. 

However, some of the PSC associations may still be 
chance alignments, and in addition the 100-/im-only sources 
could be cirrus associated with the clouds rather than em- 
bedded YSOs. So two 'control samples' were produced by 
offsetting the core positions by first 2 degrees and then 5 
degrees in declination, and selecting the MAS' point sources 
associated with the new positions. In effect this produces 
two false populations of cores, one of which is still located 
in the clouds (the 2-degree offset population), and one of 
which is located outside of the clouds. This technique was 
chosen because it was thought that it would introduce the 
least bias in the estimate of the numbers of associations that 
were purely chance alignments. 

In the control sample produced from a 2-degree dis- 
placement we found ten sources. Four were detected only at 
100/im, and five sources were detected at 12 and 25fim - 
all of which were brighter at 12^m than 25^m. These are 
thus thought to be main sequence field stars - one had been 
positively identified as a B star. The one remaining source, 
which was detected at all wavelengths, was previously cat- 
alogued as a galaxy in the Upsalla General Catalogue of 
Galaxies. In the sample produced at a displacement of 5 
degrees there were also ten sources, four of which were 100- 
^m-only sources. Six of the sources were detected at both 
12 and 25 /im and were brighter at 12 fim than 25 nm, and 
again are probably field stars. 

Hence we find no evidence for a significant difference 
between the two control samples, either due to the increased 
displacement, or to the location being inside or outside the 
clouds. Likewise, we see that the number of foreground stars 
remains roughly constant both in the control samples and 
the real sample. Of the remaining 15 PSC associations in 
the real sample, we conclude that six are probably chance 
alignments. Of the remainder, up to half may be 100-/im- 
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only cirrus sources - i.e. the point source might be the cloud 
core itself (see, for example: Benson & Myers 1989; Reach, 
Heiles & Koo 1993; Bourke et al. 1995b). Hence we estimate 
that the number of PSC sources, which are embedded YSOs 
and are associated with our sample of 60 cores, is five. 



4 COMPARISON WITH PREVIOUS CORE 
CATALOGUES 

4.1 Densities 

We can compare our new catalogue of molecular cloud cores 
with catalogues produced by previous authors using different 
methods. For example, Myers et al. (1983) surveyed 90 cores 
in C^*0 and ^'^CO. Using these observations they showed 
that the cores' C^*0 optical depths, as estimated from the 
ratio of C^*0 brightness to ^"^CO brightness, was reason- 
ably tightly distributed around a mean of 0.35 to 0.45. They 
found that the C^*0 column density, A'^is, had a mean of ~ 
1.6 X 10^^ cm"^. This led them to estimate a typical N{H2) 
column density of ~ 10^^ cm"'^. 

Clemens & Barvainis (1988) also carried out a survey of 
molecular cloud cores, and the IRAS images of these cores 
were studied by Clemens, Yun & Heyer (1991). They found 
a typical 100-/im optical depth for their sample of 2.5 x 
10~*. From this we can estimate a column density of N(H[2) 
~ 5.0 X 10^^ cm^'^, using the above equations. Wood et al. 
(1994) carried out a survey of the ISSA data for some known 
star-forming regions, and produced a catalogue of molecular 
cloud cores. They found that their cores had a mean column 
density N(H2) ~ 4.5 x 10^^ cm'^ 

The mean volume density of the cores in each of these 
surveys can also be estimated. Myers et al. (1983) found a 
typical volume density in their cores of n{H2) ~ 8 x 10^ 
cm~^. Wood et al. (1994) did not quote a typical volume 
density but a value can be estimated from their column den- 
sity if a typical radius for the cores is known. The optimum 
value to use is complicated by the fact that they defined 
cores to be areas with visual extinction A„ > 4. This leads 
to the inclusion of several very large 'cores' which we would 
define as 'clouds' (their largest 'core' is 329 pc^). This skews 
their mean to a value much larger than the median. We 
therefore take the the second quartile boundary of radius, 
which is ~ 0.5 pc in their sample. This leads to an estimate 
for the typical number density of 3 x 10^ cm""'. 

The typical number densities of the Clemens & Barvai- 
nis (1988) cores can also be estimated if the typical size of 
the cores is known. Clemens et al. (1991) claim that the typ- 
ical radius of the cores is 0.35pc. This was calculated from 
the mean solid angle of the cores and an assumed distance 
of 600 pc. However, this distance estimate is somewhat un- 
certain. It was derived originally by Clemens & Barvainis 
(1988) from two considerations: firstly that the cores were 
generally within 12 degrees of the Galactic Plane; and sec- 
ondly that the cores had LSR velocities of between and 10 
kms"'^. 

However, Bourke et al. (1995a) argued that, because 
these cores are seen in extinction against the background 
stars of the Galactic Plane, the survey is biased towards de- 
tecting cores near the Plane. Hence they derive an estimate 
of the typical distance to the cores of 300 pc. Using this 



Table 4. Summary of mean core properties of previous surveys 



of molecular cloud 


cores compared to the present work (JWT). 


See text for details. 










Survey 


N{H2) 


n(H2) 


Protostellar 


r 




Xl021 


xlO^ 


Percentage 


xlO« 




_2 

cm 


cm~^ 




yrs 


Bourke et al.(l) 


16 


30 


61±10 


0.6±0.15 


Myers et al. 


12 


8 


44±7 


1.3±0.3 


Bourke et al.(2) 


5.0 


10 


36±7 


1.7±0.4 


Clemens et al. 


5.0 


5 


27±4 


2.8±0.6 


Wood et al. 


4.5 


3 


23±3 


3.3±0.6 


JWT 


3.8 


2 


9±f 





value a typical radius of 0.175pc is derived, and hence a vol- 
ume density of n{H2) ~ 4.8 x 10^ cm~"^ is calculated for this 
sample. This is significantly higher than the value quoted by 
Clemens et al. (1991), mainly due to the different distance 
assumed, but also partly because we have followed Wood et 
al. (1994) in calculating column density from 100-/^m optical 
depth. Lemme et al. (1996) studied a subset of the Clemens 
& Barvainis (1988) cores and reached a similar conclusion: 
namely that Clemens et al. (1991) may have underestimated 
the typical density of the cores in the Clemens & Barvainis 
(1988) sample. 

Bourke et al. (1995b) presented data for two samples 
of molecular cloud cores: one sample consisted of isolated 
Bok globules, and the other was a set of cores in more ex- 
tended clouds. These had mean column densities of 5 x 10^^ 
cm~^ and 1.6 x 10^^ cm~^ respectively, and volume densi- 
ties of 10* cm"'' and 3 x 10* cm~^ respectively. Hence, it can 
be seen that each of these samples has selected cores with 
somewhat different properties. We here label the Bourke et 
al. (1995b) extended sample Bourke(l), and the Bok glob- 
ule sample Bourke(2). All of the column and volume density 
estimates have associated errors from a number of causes. 
These include chiefly the assumed fractional abundance of 
the different tracers used in each set of observations. We 
believe that the values quoted are accurate relative to one 
another to within 20 per cent. Table 4 summarises the mean 
properties of each of the core samples we have discussed. 



4.2 Protostellar content and pre-stellar life-time 

A useful parameter in the study of dense cores is the frac- 
tion of cores in a given sample that contain protostars or 
YSOs. This parameter can be used to estimate a mean sta- 
tistical life-time for the sample. This method was first used 
by Beichman et al. (1986), who studied the embedded YSOs 
within the core sample of Myers et al. (1983). They found 
that 35 cores had IRAS sources meeting the colour selection 
criteria of embedded YSOs and 43 had no embedded IRAS 
sources. 

This method was also used by Wood et al. (1994) 
with slightly different selection criteria to discard foreground 
main sequence stars. They found that 59 out of the 255 cores 
in their sample had at least one embedded source. We car- 
ried out the same test for the cores in the Clemens & Bar- 
avainis (1988) sample, and found 65 cores out of 248 have 
embedded IRAS sources (using the Beichman et al. selection 
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Figure 5. Plot of inferred pre-stellar lifetime against column den- 
sity, for each sample discussed in the text. A best- fit exponential 
of T oc e^^" Z'^) is also plotted. The reduced chi-squared error 
on this fit is 1.15. 



criteria). Bourke et al. (1995b) found that 27 out of their 76 
Bok globules had IRAS point sources satisfying the Beich- 
man criteria, while 36 out of 59 of their cores in extended 
clouds had embedded sources. 

Beichman et al. (1986) used the percentage of cores with 
embedded sources to estimate the lifetime of a core without 
an embedded YSO by comparing it with the life-time of the 
embedded YSO phase. They found the life-time of a starless 
core in this way to be a ~ lO'' years. This was based on 
assumptions relating to the time taken for a star to accrete 
its total mass, and the time for it to become visible. The 
uncertainty in this figure is probably about a factor of 2, 
but this does not affect the following statistical arguments, 
it would simply move the absolute time-scale as a whole by 
a factor of 2 (i.e. the absolute values of the vertical axes in 
Figures 5 and 6 can move up or down by a factor of 2, but 
the relative positions of the data-points do not move - see 
below) . 

Following this method, we here infer a statistical esti- 
mate of the pre-stellar lifetime, r, of the cores studied in 
each of the samples discussed above. We also make the as- 
sumption that in each survey the cores without embedded 
sources will go on to form protostars in their centres. The 
cores with embedded sources are assumed to have an aver- 
age lifetime which is the same in each sample and equal to 
the embedded YSO time-scale derived by Beichman et al. 
(1986). This is simply expressible as: 

No. of cores without embedded sources « 

r — — — — ■ — — — ■ X 10 years, 

JNo. ot cores with embedded sources 

where we have taken the lifetime of cores with embed- 
ded sources to be 10^ years as discussed above (c.f. Ward- 
Thompson et al. 1994) . r was calculated for each core sample 
discussed in section 4.1, and is listed in Table 4. 

The fraction of cores with embedded sources has ran- 
dom errors, and in the catalogues where the number of cores 
with embedded sources is large, the 1 a error is simply the 
square root of the number of cores with embedded sources. 
The number of cores in the sample presented in this paper 
has a low number of cores with embedded sources and there- 
fore the uncertainty in the pre-stellar lifetime is dominated 



Figure 6. Plot of the inferred pre-stellar lifetime against volume 
density, for each sample discussed in the text. A best-fit power- 
law of r oc n{H2)~^'^^ is plotted as a solid line. The reduced 
chi-squared error on this fit is 1.0. Functions of r oc n{H2)~^'^ 
and r cx: n(H2)~'^''^^ are also shown, as dotted and dashed lines 
respectively. 
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Figure 7. Plot of the column density against volume density, 
for each sample discussed in the text. A best-fit power-law of 
N{H2) oc n{H2)^'^^ is plotted as a solid line. 



by the systematic effects discussed in section 3.2 above. 
These errors are quoted in Table 4. 

It is apparent from Table 4 that the percentage of cores 
with embedded sources increases with both the mean col- 
umn density and volume density of the cores. Hence the pre- 
stellar core lifetime decreases with both column and volume 
densities. The results are plotted in Figs. 5 & 6. 

Figure 5 shows prestellar lifetime versus column density. 
Each of the points represents one of the data sets listed 
in Table 4 (Bl signifies Bourke et al.(l) etc.). We fitted 
a relation to the data of the form r cx exp(10'^^/A''). This 
is shown as a solid line. Figure 6 shows pre-stellar lifetime 
versus volume density (using the same labelling as in Figure 
5). Power law fits to the data points were carried out, and 
the best-fit was found to be r oc n{H2)~°-^^^°-'^^ , which is 
shown as a solid line on Figure 6. 

Figure 7 plots column density against volume density 
for each of the core data sets in Table 4, using the same 
labelling once again. The solid line is a best fit to the data. 
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which is N{H2) oc n{H2)°-'^''"^°-^^ . This is of a similar form 
to the empirical observations usually referred to as Larson's 
relations (Larson 1981), although with a somewhat differ- 
ent slope. If we make the assumption that each catalogue 
is a representitive subset of all star-forming clouds, and the 
different membership of each subset is simply due to the 
characteristic density sensitivity of the different tracers used 
in each study, then we can compare these empirical results 
with theoretical predictions. 

4.3 Comparison with theory 

We can compare the dependence of prestellar time-scale on 
volume density seen here, with the ambipolar diffusion time- 
scale. However, it must be noted that the densities plotted 
in Figure 6 are the mean densities for each core sample, 
volume-averaged over the whole core in each case. They are 
not the central densities of each core. Hence care must be 
taken when comparing Figure 6 with ambipolar diffusion 
models, which usually plot central density versus time-scale. 

The ambipolar diffusion time-scale is proportional to 
the ionisation fraction, Xi- The ionization of the gas can be 
caused both by ultra violet radiation and cosmic ray ioni- 
sation. The canonical form for the relation between cosmic 
ray ionisation and volume density is usually taken to be a 
power-law. For example, Mouschovias (1991) uses: 

X^ cx n(^/2)-0•^ 

which leads to the derivation of the ambipolar diffusion time- 
scale relative to density of: 

TAD OC n(H2) 

This behavior is plotted as a dotted line in Fig. 6. When 
additional factors are included (such as chemistry, multiple 
charge carriers etc.), the volume density has a slightly differ- 
ent influence on recombination rates, and hence ionisation 
levels. For example, McKee (1989) suggests: 

/ rr \-0.75 

Xi OC n{H2) 
leading to: 

TAD 0Cn(ii-2)~''-''. 

This is plotted as a dashed line on Figure 6. It can be seen 
that the steeper relation is a closer match to the data, and 
is in fact consistent to within 1 a. 

The fact that the lowest density point on Figure 6 lies 
above the best-fit line suggests that at low densities an ad- 
ditional effect may be starting to become significant. This is 
somewhat tentative, but does have a possible theoretical ex- 
planation. For example, McKee (1989) treats ionization and 
recombination in some detail, and incorporates the role of 
metals, and UV penetration. He derives an expression for the 
star formation time-scale which is dependent on both vol- 
ume density and column density (see his equation 4.5, and 
his fig. 1). He also shows that UV ionisation affects the star 
formation timescale, yielding an exponential dependence on 
column density (McKee equation 5.2) of the form: 

(Na/N) 

TSF QC ', 

where tsf is the star formation timescale and Nc is the 
critical column density, tsf is not the same as the timescale 



we have plotted in Figure 5, since McKee was referring to the 
timescale in which an entire molecular cloud is converted to 
stars and we are considering the timescale for the dense cores 
within clouds to form stars. Hence his timescales are roughly 
an order of magnitude greater than ours. However, we may 
be observing a similar exponential behaviour. The constant 
in the exponent, Nc, was derived by McKee to be 1.6 x 
10^^ cm~^ (deduced from his critical extinction estimate of 
Av=16), which is consistent, to within the errors, with the 
value of 1.0 x 10^^ cm~^ that we derive in Figure 5. 

The three parameters of column density, volume den- 
sity and lifetime are all inter-dependent. Hence Figures 5 & 
6 are not independent. Therefore the empirical fits which 
we have applied to these plots are not strictly separable in 
the manner we have adopted. However, we used this ap- 
proach for the sake of clarity and of demonstrating the po- 
tential underlying physical processes. From these data we 
cannot make definitive statements about ionisation mecha- 
nisms, but we have shown that the data, not only from our 
current survey, but also from those earlier surveys that we 
have discussed above, are all consistent with a picture in 
which the ionisation levels in molecular cloud cores regulate 
the star formation timescale. 



5 CONCLUSIONS 

We have presented a catalogue of 60 cores situated in 
medium opacity molecular clouds, with the aim of broaden- 
ing the range of physical environments in which star forma- 
tion has been studied. The catalogued cores typically have 
lower column densities and volume densities than previously 
studied samples and a lower fraction of the cores have formed 
stars. We found a clear trend for cloud cores to form stars 
more rapidly with increasing volume and column density. We 
hypothesised that this can be interpreted in the framework 
of ionisation-regulated star formation. 
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